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ABSTRACT 

At modest radii from the centre of galaxy clusters, individual galaxies may be in- 
falling to the cluster for the first time, or have already visited the cluster core and are 
coming back out again. This latter population of galaxies is known as the backsplash 
population. Differentiating them from the infalling population presents an interesting 
challenge for observational studies of galaxy evolution. To attempt to do this, we as- 
semble a sample of 14 redshift- and spatially-isolated galaxy clusters from the Sloan 
Digital Sky Survey. We clean this sample of cluster-cluster mergers to ensure that the 
galaxies contained within them are (to an approximation) only backsplashing from the 
centre of their parent clusters and are not being processed in sub-clumps. By stacking 
them together to form a composite cluster, we find evidence for both categories of 
galaxies at intermediate radii from the cluster centre. Application of mixture mod- 
elling to this sample then serves to differentiate the infalling galaxies (which we model 
on galaxies from the cluster outskirts) from the backsplash ones (which we model on 
galaxies in the high density core with low velocity offsets from the cluster mean) . We 
find that the fraction of galaxies with populations similar to the low velocity cluster 
core galaxies is / = — Q.052R/ R V i r i a i + 0.612 ± 0.06 which we interpret as being the 
backsplash population fraction at 1 < R/Rviriai < 2. Although some interlopers may 
be affecting our results, the results are demonstrated to be in concordance with earlier 
studies in this area that support density-related mechanisms as being the prime factor 
in determining the star formation rate of a galaxy. 
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1 INTRODUCTION 

During the course of its lifetime, a galaxy may travel through 
several different environments and be affected by different 
physical mechanisms in each of them. As a very simplistic, 
illustrative example, consider an 'average' test galaxy that 
gets accreted on to a filament of galaxies. As the galaxy flows 
along the filament, it may undergo some pre-processing due 
to first time harassment (Moore et al. 1996) caused by the 
increase in local galaxy density (Porter et al. 2008). The fil- 
ament will eventually funnel the galaxy to anisotropically 
fuel the growth of a cluster of galaxies (Pimbblet 2005; see 
also White, Cohn & Smit 2010). If the star formation rate 
of the galaxy has not already undergone significant trans- 
formation due to (slow) gas starvation (cf. van der Wei et 
al. 2010; Balogh, Navarro, & Morris 2000; Larson, Tinsley, 
& Caldwell 1980), then the cluster environment will almost 
certainly act to truncate it to much lower levels (e.g. Lewis et 
al. 2002; Gomez et al. 2003; Kauffman et al. 2004; Pimbblet 
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et al. 2006) as higher pressure and higher density mecha- 
nisms such as ram pressure stripping (Gunn & Gott 1972) 
and tidal interaction with the cluster gravitational poten- 
tial (Byrd & Valtonen 1990; Rose et al. 2001) become more 
important. Additionally, ram pressure stripping may also in- 
duce a final star star formation episode in a galaxy to remove 
any remaining gas (Smith et al. 2010a). 

These mechanisms may also change the morphology of 
the galaxy by (e.g.) perturbation of spirial arms (Quilis, 
Moore, & Bower 2000) and the colour will get progressively 
redder as the star-formation rate declines and the bluer stel- 
lar population dies off. Such processes are generally posited 
to contribute to the well-known morphology- density and 
density-star formation rate relationships (e.g. Dressier 1980; 
Smith et al. 2005; Holden et al. 2007). Taken together, the 
origin of these relations appear dependant on a variety of 
physical mechanisms that act on distinct galaxy masses, at 
different epochs, and over certain time-scales. 

Attempting to analyze the star formation rates of such 
galaxies in order to extract more general conclusions about 
the driving mechanisms behind galaxy evolution and the 
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morphology-density relation is not without its trials. For 
example, the work of Porter et al. (2008) examines the star- 
formation rate of galaxies as a function of cluster-centric ra- 
dius along the vector of a filament going away from a cluster 
core. They suggest there is a peak in star-formation rate at a 
few Mpc from the cluster core along filaments driven mostly 
by light from dwarf galaxies - this is something that would 
not be seen if one simply averages all vectors away from the 
cluster core, as works such as Lewis et al. (2002), Gomez et 
al. 2003, and Pimbblet et al. (2006) do. 

One further complication to addressing the changes in 
star formation rate (whether radially averaged, or chosen 
along the vector of a filament as in Porter et al. 2008) is 
the presence of backsplash galaxies. Backsplash galaxies are 
those galaxies that have already entered the core regions 
of a galaxy cluster and have subsequently come back out 
again (Gill, Knebe, & Gibson 2005; see also Ludlow et al. 
2009; Warnick et al. 2008; Aubert & Pichon 2007). The exis- 
tence of backsplash galaxies is not unexpected: Balogh et al. 
(2000) find that 50±20% of particles in the region of R 200 
to 2i?2oo from a cluster have passed within 7?2oo- Mamon 
et al. (2004) further compute that the maximum distance 
that a backsplash particle has is 2.5i?ioo, and for a galaxy is 
1.7Jiioo- This is supported by Moore et al. (2004) who com- 
pute that half of all halos between R v i r iai and 2R v i ria i have 
passed through their parent cluster, many directly through 
the high density core regions of clusters Interestingly, Knebe 
et al. (2008) show that the backsplash population is pre- 
dicted to exist in more than one cosmology. 

The implication of these findings is observationally con- 
cerning: to understand the origin of (e.g.) the morphology- 
density relation and the mechanisms that drive galaxy evolu- 
tion, we should in principle objectively disentangle the back- 
splash galaxy population from those that are infalling to the 
cluster for the first time. 

Gill et al. (2005) provide a succinct observational 
method to detect the presence of a backsplash population in 
galaxy clusters. Given that the backsplash population will 
have a distinct velocity distribution, observers need simply 
to examine the distribution in a suitably large radius interval 
from the cluster centre, say between 1 and 2 virial radii. A 
backsplash population will exhibit a centrally peaked distri- 
bution for this interval, whereas a pure infalling population 
will possess a higher modal value. 

A number of observational studies have now confirmed 
the existence of a backsplash population in galaxy clusters 
(e.g. Aguerri & Sanchez-Janssen 2010; Smith et al. 2010b; 
Sato & Martin 2006; Pimbblet et al. 2006; Rines et al. 2005; 
Sanchis et al. 2004). Yet, although these studies strongly 
suggest that the backsplash population exists, they have not 
in detail disentangled the backsplash population from the 
infalling one, although using galaxies with different spectral 
features appears to hold much promise for this (see Pimbblet 
et al. 2006; Rines et al. 2005). 

This work revisits the backsplash concept and (observa- 
tionally) attacks the question of how the backsplash fraction 
is expected to vary with radius from the cluster centre. In 
Section 2, we construct a sample of isolated galaxy clus- 
ters from the Sloan Digital Sky Survey to work with. This 
sample should be free of significant nearby structures and 
sub-clumps that may be acting to pre-process before they 
reach the centre of their parent galaxy clusters. In Section 3 



we observe the presence of the backsplash population in our 
sample. We apply a mixture model to deduce how the frac- 
tion of backsplash galaxies varies with radius from the centre 
of clusters and compare our result with previous studies in 
this area. In Section 4 we generalize our results and discuss 
a number of caveats about our work before summarizing the 
main findings in Section 5. 

Throughout this work, we adopt Ho = 100 
/iiookms -1 Mpc -1 , Q. M = 0.27, and Q. va cuu m = 0.73. 



2 SAMPLE 

For this work, we wish to use a sample of galaxy clusters that 
are free from significant subclustering that may skew some 
of the earlier observational results. For instance, Pimbblet 
et al. (2006) make use of very X-ray bright, massive clusters, 
some of which exhibit significant substructure (e.g. A1664; 
A3888). The issue is that inside these sub-clusters, galax- 
ies may have been already effected by various evolutionary 
mechanisms as if they had fallen into a cluster already. We 
wish to avoid this scenario as much as possible. Fortuitously, 
such an ensemble already exists: Plionis et al. (2009) have as- 
sembled a large sample of galaxy clusters from Abell, Corwin 
& Olowin (1989) that are 'clean' of merging and interacting 
clusters. This sample should be almost ideal for studying a 
backsplash population with, since in principle it should be 
free from coherent infall (e.g. large groups and sub-clusters 
of galaxies) and relatively isolated in space. 

From the sample of Plionis et al. (2009) we select only 
those clusters that are contained within Sloan Digital Sky 
Survey (SDSS; Abazajian et al. 2009) boundaries and ap- 
pear in the SDSS C4 catalogue (Miller et al. 2005) to study 
further - this is a matter of convenience given SDSS has 
good photometry down to a reasonably faint limit for nearby 
galaxy clusters, coupled with quality spectra. Moreover, by 
using SDSS as our exclusive data source, we try to actively 
avoid any biases that might be present in ensembles from 
inhomogeneous data sources. These clusters are detailed in 
Table 1. Given that all of them have been detected by the 
C4 cluster finding algorithm (Miller et al. 2005), they are 
not chance superpositions along the line of sight as a non- 
negligible fraction of Abell clusters have been shown to be 
(cf. Pimbblet et al. 2010; Frenk et al. 1990; Lucey 1983). 

From this sample, we discard A1452. The reason for 
doing so is that in the C4 catalogue of Miller et al. (2005), 
this cluster is noted to have a high probability of having sub- 
clustering, as determined from a Dressier & Shectman (1988) 
test, despite its inclusion in the Plionis et al. (2009) sample. 
Additionally, A0095 is also flagged as not being highly iso- 
lated in redshift space (based on Miller et al. 2005's 'struc- 
ture contamination flag'). We therefore discard this cluster 
as well. This leaves us with a sample of 14 galaxy clusters 
to work with. 

We detail the global parameters of these clusters in Ta- 
ble 1 (e.g. mean redshift, ~cz\ velocity dispersion, a cz ; X-Ray 
luminosity in the 0.1-2.4 keV band, L x )- We note that the 
values presented in Table 1 (i.e. sourced from Miller et al. 
2005) compare very favourably to other, sometimes more 
extensive and deeper, studies of galaxy clusters (e.g. A1650 
also appears in the Pimbblet et al. 2006 sample and exhibits 
completely consistent values for a cz as used here) . From the 
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values of a cz , we compute the virial radius (R v i r iai) for each 
cluster from the expression given by Girardi et al. (1998): 
-R200 ~ Rviriai = 0.002<Tcz (where -R200 is the radius at which 
the mean interior density is 200 times the critical density 
and is well approximated by Rviriai)- Although other ex- 
pressions exist for similar radii (e.g. Carlberg et al. 1997 give 
Rviriai ~ #200 = 3 1/2 a cz /WH(z), assuming an isothermal 
sphere), this expression will suffice for ultimately placing our 
clusters on to a common metric. 

We note that the sample spans a factor of 1000 in X-Ray 
luminosity; from clusters that are often thought of as being 
amongst some of the most massive in the local Universe 
(A1650; cf. Pimbblet et al. 2006) to those whose luminosity 
is comparable to rich groups, or lower. The velocity disper- 
sion of these clusters, however, is more in-line with that of 
intermediate to massive clusters (see Table 2 of Pimbblet et 
al. 2006) with a median value of 849 kms" 1 (see also Ortiz- 
Gil et al. 2004 who discuss the peculiarities of the Lx-<J CZ 
relationship in depth). This implies that the range in total 
mass of our clusters is likely to be little more than a factor of 
~ 2 assuming o cz is a reasonable proxy for mass. The range 
of o cz spans a factor of ~2.4. This is shortened to only a 
factor of ~ 1.6 when the highest and lowest values in the 
sample are ignored. In terms of look-back time, the clusters 
span ~1 Gyr range at maximum (from A0119 to A0655), 
or an inter-quartile range of ~0.3 Gyr - short enough that 
the time-scales for environmental processes driving galaxy 
evolution will not greatly affect the sample (e.g. Kodama & 
Smail 2001). Or put another way: the clusters should all be 
at a very similar stage in their evolution. 

Consequentially, we suggest that the sample contains 
reasonably self-similar clusters that should have comparable 
backsplash populations for the purposes of this study. 



3 BACKSPLASH IN THE COMPOSITE 
CLUSTER 

To proceed, we now extract all galaxies from SDSS that are 
within SRviriai of their cluster centres, as given in Table 1. 
The homogeneity of the sample is then exploited in order 
to construct a stacked, composite cluster. This is done by 
placing the cluster on to a standard fixed metric through 
the use of the computed Rviriai values (Table 1) and scal- 
ing each cluster by its velocity dispersion. The result of this 
stacking process is illustrated in Fig. 1 where all galaxies 
within SRviriai are plotted against their scaled velocity dif- 
ference from cz. To better define the extent of the cluster, 
we follow Pimbblet et al. (2006; see also Lerchster et al. 
2010) and utilize the cluster mass model of Carlberg, Yee 
& Ellingson (1997; see also Carlberg et al. 1997; Rines & 
Diaferio 2010 and references therein). Only galaxies inside 
the mass model's caustics (Fig. 1) are considered members 
and used in our subsequent analysis. The relative isolation of 
the composite cluster in both a radial and recession velocity 
sense can clearly be seen from Fig 1, even out to modestly 
large radii from the cluster centre; thereby validating our 
approach of using a reduced Plionis et al. (2009) sample. 

Next, we split the sample in to a number of radial 
bins. Following Gill et al. (2005), we create histograms of 
\Acz\/a cz in each of these radial bins and display these dis- 
tributions in Fig. 2. A pure backsplash population should 
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Figure 1. The composite, stacked cluster plotted as the differ- 
ence in velocity normalized to each cluster's individual velocity 
dispersion as a function of virial radius. The solid curves denote 
the 3<r contour (caustics) of the mass model of Carlberg, Yee & 
Ellingson (1997). Only galaxies inside these caustics are consid- 
ered members and used in the subsequent analysis. 



appear centrally peaked in these plots, as per the predictions 
of Gill et al. (2005), whilst an infalling population will be 
peaked in the 0.3 < \Acz\/a cz < 0.5 region (which we note is 
contained inside the mass model limits at all radii explored; 
Fig. 1). Indeed, the infalling population can be expected to 
have approximately a factor of two greater relative velocity 
than the backsplash ones (Gill et al. 2005) which makes the 
two populations kinematically distinct. Gill et al. (2005) re- 
port that if all galaxies at modest radii (~ few Rviriai) are 
infalling, then their absolute velocity has a modal value of 
\cz — cz\ w 400 kms -1 . For clusters such as those used in 
this work (and also Pimbblet et al. 2006 & Rines et al. 2005) 
this corresponds to approximately 0.3 < \Acz\/a cz < 0.5. 
Fig. 2 shows evidence for both types of populations being 
present, similar to previous observational studies (Rines et 
al. 2005; Pimbblet et al. 2006). At the centre of the compos- 
ite cluster (R/ Rviriai < 0.5), the distribution is markedly 
centrally concentrated in \Acz\/a cz terms. The apparent de- 
gree of central concentration decreases with increasing ra- 
dius. In the zone where we expect a mixed population (e.g., 
1.0 < R/ Rviriai < 1-5), we see that the mode of the distribu- 
tion sits in the range 0.3 < \Acz\/a cz < 0.5, as would be ex- 
pected from an infalling population (Gill et al. 2005), along 
with many galaxies under this range (i.e. the backsplash 
population). As we go out to higher radii, the distinction of 
the modal value lessens. Furthermore, at large radii, we do 
not expect any backsplash galaxies to be present (Mamon 
et al. 2004). 
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Table 1. The sample of Abell clusters used in this work. The mean redshift (cz), and velocity dispersion (cr cz ) are sourced from Miller 
et al. (2005). The coordinates specify the Vizier position of the cluster. The X-Ray luminosities (Lj) are sourced, if available, from the 
Base dc Donnccs Amas de Galaxies X (BAX; webast.ast.obs-mip.fr/bax; Sadat ct al. 2004) which in turn derives its data from other 
literature sources (Popcsso ct al. 2007; Rciprich & Bohringcr 2002; see also Ebeling et al. 1996). The virial radius {R v i r i a i) is computed 
from cr cz \ see text for details. 
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^A0095 is removed from our sample on the basis of a high structure contamination 
flag (Miller ct al. 2005). 

tA1452 is removed from our sample since Miller et al. (2005) indicate a high probability 
of subclustcring. 



3.1 Specific Star Formation Rates 

In addition to the difference in the dynamics of the infall 
and backsplash population, a difference in star formation 
histories may also be present. If the backsplash population 
has already experienced the hostile high density core regions 
of a galaxy clusters (cf. Moore et al. 2004), then it fol- 
lows that their star formation rate distribution should reflect 
this. This provides us with a method to potentially deduce 
the minimum fraction of backsplash galaxies present in the 
mixed population (Fig. 2). The reason that it will be a min- 
imum fraction is that not 100% of the backsplash galaxies 
will have experienced star-formation truncation compara- 
ble to core-region galaxies - some may have more fortunate 
orbits around the cluster centre and retain an appreciable 
star formation rate. Indeed, this is supported by Rines et 
al. (2005) who demonstrate that some backsplash galaxies 
must be emission line galaxies. Hence the status of a given 
galaxy as being a backsplash one may not be the prime de- 
terminant of their star formation. Therefore in the analysis 
below we examine the specific star formation rate distribu- 
tion of a large number of galaxies at various locations inside 
the Carlberg, Yee & Ellingson (1997; Fig. 1) cluster mem- 
bership caustic. 

If the galaxy population located at 0.3 < \Acz\/a C z < 
0.5 is considered to be a mixture of both infalling and back- 
splash galaxies, then we should be able to model it as a 
superposition of both and determine the how the fraction of 
backsplash galaxies varies as a function of radius from the 
cluster centre (cf. Rines et al. 2005; see also Pimbblet et al. 
2006). 

Arguably the simplest population to get a handle on is 
the infalling population. The simulations of Mamon et al. 



(2004) demonstrate almost no backsplash galaxy will reside 
at > 1.7rioo- We therefore assume that at radii greater than 
2R v iriai (Fig. 2), the galaxy population consists of approxi- 
mately pure first-time infallers. For galaxies in the range 2-3 
Rviriai, we find the logarithm of the mean specific star for- 
mation rate (SSFR; in units of MQyr^ 1 M*^ 1 ; see Brinch- 
mann et al. 2004 and www.mpa-garching.mpg.de/SDSS for 
details of this parameter's derivation) is —10.99 ± 0.88 with 
an inter-quartile range of 1.79. If we restrict this sample to 
0.5 < \Acz\/a C z < 2.0 (i.e. well outside the centrally con- 
centrated peak one might expect for backsplash galaxies, 
and outside the plausible maximum backsplash radii) the 
result does not vary considerably: —10.97 ± 0.92 with an 
inter-quartile range of 1.78. 

We next turn to the core population. There may be 
some ambiguity in defining the core population due to infall 
interlopers at R < R viria i whose three dimensional radius 
from the cluster centre is > R v i r iai (Diaferio et al. 2001; 
Rines et al. 2005). The fraction of such interlopers is ex- 
pected to be low, however (perhaps 10 to 20%). Therefore 
in order to get a handle on the specific star formation rates 
in the cluster core, we restrict ourselves to not only small 
radii from the cluster centres (R < 0.5R viria i), but also to 
those galaxies that are highly centrally peaked (explicitly: 
\Acz\/a cz < 0.3). This is reasonable given the predictions 
of the large number of backsplash galaxies and particles 
that have visited this zone (Balogh et al. 2000; Moore et 
al. 2004). For these galaxies, the mean SSFR is found to be 
logio(SSFR) = —11.50 ± 0.73 with an inter-quartile range 
of 0.67. 

In Fig. 3, we demonstrate how the specific star forma- 
tion rate of the mixed population (0.3 < \Acz\/a cz < 0.5) 
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Figure 2. Histograms of |Acz|/o- cz in the composite cluster, split 
in to radial bins. The mode of the distribution is expected to 
reside between the vertical dotted lines for a pure population of 
first time infallers. A mixture of both backsplash and infalling 
populations can be seen at low to intermediate radii. 

varies with distance from the cluster centre. In line with 
other results (cf. Gomez et al. 2003), we observe a decline of 
SSFR at the centre of the composite cluster. Importantly, 
from Fig. 3, it appears to be the case that the variation of 
the mixed galaxy population with radius can be effectively 
thought of as a superposition of the core and infalling pop- 
ulation. 

3.2 Mixture Modelling 

We have thus far not attacked the central question of the 
relative fraction of backsplash galaxies with radius in the 
composite cluster. To do this, we will need to evaluate what 
relative proportions of galaxies in the composite cluster have 
SSFRs that are consistent with either the core population 
(i.e. consistent with having their star formation rates trun- 
cated due to interactions with the hostile cluster core) or the 
infalling population as computed above (cf. Fig 3). This can 
be readily achieved using a mathematical mixture model. 
One of the most frequently used mixture models in galaxy 
cluster analysis is Kaye's Mixture Model (KMM; Ashman 
et al. 1994). Although KMM is extensively used in the lit- 
erature to segregate substructure in clusters (e.g. Pimbblet 
et al. 2010; Owers et al. 2009), it is very much a general 
purpose tool that can be applied to the present situation. 
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Figure 3. Variation of the mean and inter-quartile SSFR range of 
galaxies in the range 0.3 < \Acz\/a cz < 0.5 with radius from the 
cluster centre (solid lines). The inter-quartile range for first time 
infallers (R > 2i?„; r ; (1 ;) and the core population (R. < 0.5R v i r i a i 
and \Acz\/cr cz < 0.3) are displayed for reference as the horizontal 
dashed and dotted lines respectively. We can model the change 
in radius as a mixture of these infalling and core populations. 



As Gaussian inputs to the KMM algorithm we supply 
the means and standard deviations of the SSFR for the core 
and infalling population that we found above. Then for ra- 
dial bins in the 0.3 < \Acz\/a cz < 0.5 slice, we implement 
KMM to partition the data and inform us of the relative 
fractions of galaxies that are consistent with either the low 
velocity offset, core region population or infalling population 
(as defined above). The result of this analysis is presented 
in Fig. 4. Although there is scatter in the points, the trend 
shows a downward trajectory with radius. In the mixed pop- 
ulation radial regime (1 < R/R V i r iai < 2), we find that 
the fraction of galaxies with SSFRs similar to the low ve- 
locity offset cluster core region that we are interpreting as 
backsplash galaxies varies from 0.68 ± 0.09 (at R V i r iai) to 
0.55 ±0.09 (at 2R vlrlal ). 

We note that the relationship presented in Fig. 4 tells 
us only about the fraction of galaxies that have SSFRs con- 
sistent with the low velocity core galaxies. If we interpret 
these as backsplash galaxies (i.e. by assuming they have vis- 
ited the hostile cluster core and had their rates truncated), 
we caution that the relation should only be applied at radii 
where mixed backsplash and infall galaxy populations are 
expected: especially 1 < R/R V iriai < 2. Outside of this 
regime (>> 2R/R viria i), there are only low (or zero) back- 
splash galaxies expected (cf. Mamon et al. 2004) whereas 
below R V iriai , we are only detecting the core population that 
we modelled the backsplash galaxies on in the first instance. 

We now ask whether our calculated backsplash values 
agree with previous works? One of the few predictions on 
this value comes from Gill et al. (2005; their Fig. 8). That 
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Figure 4. Radial fraction of galaxies in the 0.3 < | Acz\/a C z < 0.5 
slice consistent with the distribution of SSFR of the low velocity 
offset core galaxies (Fig. 3) that we interpret as the backsplash 
fraction. Each point contains 50 galaxies and the errorbar is the 
poisson error for this. The solid line shows the linear line of best fit 
(/ = — 0.135-R/-R„i r i a ; + 0.823) to the data points together with 
Iff errors on the fit's intercept (dotted lines; ±0.086). The ver- 
tical dot-dash lines denote the expected mixed population radial 
regime (1 < R/R v i r ial < 2) where these fractions are expected to 
be of maximum utility (i.e. at >> 2R/ R viria i there are low, to 
zero backsplash galaxies expected) . These fractions are consistent 
with previous estimates in other works. 

prediction (reproduced in Rines et al. 2005; their Fig. 8, top 
panel) show that the backsplash fraction is expected to be 
« 0.62 at 1.4 < R/r 2 oo < 2.8 and \Acz\/a 2 oo ~ 0.4. This is 
in agreement with our best fit line (Fig. 4). 



4 DISCUSSION 

To extend the utility of our results, we show in Fig. 5 how 
the fraction varies as a function of radius for the entire clus- 
ter population - i.e. with no cut made on \Acz\/a cz other 
than they fall within the cluster membership caustic (Fig. 1). 
The form of the slope (/ = -0.052^/^™; +0.612) is much 
shallower than for the 0.3 < \Acz\/a cz < 0.5 slice (Fig. 4). 
At R v iriai, Fig. 5 suggests a fraction of 0.59 ± 0.06 galaxies 
are backsplash, dropping to 0.52 ± 0.06 at 2Rvirial- We sus- 
pect the reason for this shallower relationship is the increase 
in the infalling interloper fraction due to making no cut in 
\Acz\/cr cz . Indeed, Rines et al. (2005) demonstrate that up 
to 60% of emission line galaxies below r 2 oo (~ Rviriai) are 
probably interlopers that have a three dimensional radius 
from the cluster of > R V i T iai- Given that such emission line 
galaxies have significantly higher values of \Acz\/a cz in clus- 
ters (Pimbblet et al. 2006; Rines et al. 2005; see also Biviano 
& Katgert 2004), we regard it as very likely that this is the 
prime cause of the change between Fig. 4 and Fig. 5. 
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Figure 5. As for Fig. 4, but with no cut made on |Acz|/ct cz 
and each point containing 500 galaxies. The observed trend (/ = 
—0.052R/ R v i ria [ +0.612) is shallower than that depicted in Fig. 4. 

From the outset, we caution that the sample of galaxies 
that we have used has not been volume limited. Although 
we may have, in principle, circumvented this by the use of 
specific star formation rates (Fig. 3), it remains the case 
that some environmental galaxy evolution mechanisms may 
preferentially operate on lower mass galaxies at different 
radii and densities than more massive ones (e.g. Penny et 
al. 2010 and references therein). However, if anything, we 
suspect that the star formation suppression that occurs in 
lower mass galaxies may be more severe than in high mass 
ones (e.g. Boselli et al. 2008), leading to a higher fraction 
of backsplash galaxies in dwarf regimes at comparable radii 
than for their more massive cousins. To test this and create 
a volume-limited sample, we must necessarily truncate the 
clusters used to the same portion of the luminosity func- 
tion. For SDSS spectroscopy of the main galaxy sample, we 
note that the limiting magnitude for ~ 100% completeness is 
r = 17.77 (Abazajian et al. 2009; Strauss et al. 2002). We use 
a simple stellar population coupled with a solar metallicity 
and 2 = 5 formation time as input to Bruzual & Chariot's 
(2003) code to predict how this magnitude limit varies as 
a function of redshift. Between our highest redshift cluster 
(A0655) and lowest (A0119), the r = 17.77 limit corresponds 
to a change in absolute magnitude of AMr ~ 2.5. Such a 
large change would mean that our lowest redshift clusters 
would contribute only very small numbers of galaxies (~l's 
to 10's of galaxies) to the composite, and overall we would 
be reduced to a small (< 100 galaxies) sample size. To alle- 
viate this issue, we consider clusters in a narrower redshift 
slice: from A1620 {z = 0.084) to A1767 (z = 0.070). This re- 
duced sample totals 6 clusters and spans a change of absolute 
magnitude of AM R ~ 0.4. For A1620 (the highest redshift 
cluster in the reduced range), the SDSS apparent magnitude 
completeness limit corresponds to an ~ M* + 2 galaxy. We 
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Figure 6. As for Fig. 4, but for the reduced, volume lim- 
ited sample of 6 clusters with no cut made on \ Acz\/a cz . Each 
point contains 50 galaxies. The fitted trend is of the form / = 
-0.097R/R virial +0.765. 

now truncate the galaxies in the other 5 clusters in the nar- 
rower redshift regime to the same limit and reproduce the 
above analysis on the reduced sample. Fig. 6 displays these 
results. At Rviriai, the backsplash fraction from the fitted 
line is 0.67 ± 0.10, decreasing to 0.57 ± 0.10 at 2R virial . We 
therefore regard these volume limited results as compatible 
(i.e. within la) with the previous results (Fig. 5). 

Lastly, we note here that changing the exact choice of 
the cluster centre to (e.g.) the brightest cluster member has 
no significant effect on the fits given in Figs. 4 through 6. 



5 SUMMARY 

The main findings of this work are summarized below. 

(1) We have detected a mixed population of infalling 
and backsplash galaxies in isolated (in both redshift space 
and the plane of the sky) galaxy clusters that are free from 
recent cluster-cluster merger activity. 

(2) We have modelled the backsplash population on low 
velocity offset core population members under the assump- 
tion that backsplash galaxies will have had their SSFRs re- 
duced from interaction with the hostile cluster core. This 
is reasonable given the predictions of the sheer numbers of 
backsplash galaxies that have visited this zone (cf. Moore 
et al. 2004). By using a mixture model, we suggest that 
this population follows the relation / = — 0.052R/ R v i r i a i + 
0.612 ± 0.06 for all cluster members. 

We explicitly note that this equation should only be 
applied where there is a mixed population expected - i.e. not 
in the very high density core, and not at the outskirts of the 
cluster where backsplash galaxies are unexpected to reside. 
Ideal radii for application is over the range ~ Rviriai to 



~ 2R V i r iai. We caution that this relation makes no account 
of infall interlopers. 

(3) Our results are in broad agreement with previous 
observational works in this area (e.g. Pimbblet et al. 2006; 
Rines et al. 2005) who do not use such idealized isolated 
samples. 

Taken together, our results support the viewpoint that 
mechanisms tied to local galaxy density should be more im- 
portant drivers of star formation than distance to galaxy 
clusters. This is in agreement with, and adds weight to, our 
earlier results (Pimbblet et al. 2006) and is opposite to (e.g.) 
Whitmore, Gilmore & Jones (1993). 

In the near future, it will be interesting to closely ex- 
amine the results of the Galaxy and Mass Assembly survey 
(GAMA; Driver et al. 2010) to see what is happening in the 
more dwarfier regimes, and the on-going investigations of 
Owers et al. (2009) who have obtained very high spectro- 
scopic completeness in a batch of clusters. 
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